HD-TVP-94-6 



QUARK-GLUON-PLASMA 
FORMATION 
AT SPS ENERGIES?i 

K. WERNERgfl 

Institut fiir Theoretische Physik, Universitat Heidelberg 
Philosophenweg 19, 69120 Heidelberg, Germany 

and 

Max-Planck-Institut fiir Kernphysik, Heidelberg, Germany 



By colliding ultrarelativistic ions, one achieves presently energy densities close 
to the critical value, concerning the formation of a quark-gluon-plasma. This 
indicates the importance of fluctuations and the necessity to go beyond the 
investigation of average events. Therefore, we introduce a percolation approach 
to model the final stage (r > 1 fm/c) of ion- ion collisions, the initial stage 
being treated by well-established methods, based on strings and Pomerons. The 
percolation approach amounts to finding high density domains, and treating 
them as quark-matter droplets. In this way, we have a realistic, microscopic, 
and Monte— Carlo based model which allows for the formation of quark 
matter. We find that even at SPS energies large quark-matter droplets are 
formed - at a low rate though. In other words: large quark-matter droplets are 
formed due to geometrical fluctuation, but not in the average event. 
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Based on extrapolations of present knowl- 
edge, one expects that at future ion-ion col- 
liders (RHIC, LHC) very high energy den- 
sities will be achieved, high enough to sup- 
port the formation of a new state of mat- 
ter consisting of deconfmed quarks and glu- 
ons. At present days accelerators for ultra- 
relativistic ions (SPS, AGS), operating at 
lower energies, the situation is less clear, 
because one reaches energy densities just 
around the critical value. In such a situ- 
ation, fluctuations become crucial, and it is 
no longer useful to discuss "average events" . 
Therefore, we investigate the formation of 
quark-gluon-matter droplets due to geomet- 
rical fluctuations. This amounts to find- 
ing domains of high energy density - higher 
than the average. For this purpose, we in- 
troduce a percolation approach for the final 
stage, r > 1 fm/c, of nucleus-nucleus colli- 
sions, the initial stage being treated by the 
VENUS string model 0. 

Straightforward extrapolations of models 
for (soft) hadron-hadron scattering turned 
out to be quite successful for hadron-nucleus 
as well as for nucleus-nucleus scattering. 
This is, however, only true for very basic 
observables like the pion multiplicity or the 
transverse energy. The investigation of rare 
processes (like A production) showed the 
limitations of the "straightforward extrap- 
olations" and the necessity to go beyond. 
The percolation approach to be discussed 
later is such a step beyond, other mecha- 
nisms have been proposed like a hadronic 
cascade |2], |3|, || or string fusion |J [7J]. 
The percolation appraoch is the first real- 
istic, microscopic, and Monte-Carlo based 
model to allow for the formation of quark 
matter. 

Let us first discuss the basic ideas of the 
percolation approach in a schematic way, 
the details and, in particular, the appro- 
priate relativistic formulation will be given 
later. Consider a snapshot at some fixed 
time r. This time shold be large enough 



Figure 1: Overlapping volumes, 
representing high energy densi- 
ties. 

(> 1 fm/c), so that the "initial stage" in- 
teractions occuring in the nucleus-nucleus 
overlap zone are finished. At the given 
time r, we consider the locations (in R 3 ) 
of all particles produced in the initial stage. 
There are, by chance, regions with high (en- 
ergy) density and such with low density. To 
be quantitative, we introduce a "critical en- 
ergy density" Eq, and look for domains with 
s > 6q. This can be done for example by 
using a grid and checking the density per 
cell. The high density domains, the con- 
nected regions with e > Eq, are referred 
to as quark matter clusters (or droplets). 
Once they have been formed, these clus- 
ters are treated macroscopically, character- 
ized by a distribution function for energy 
density Ei(x, r), momentum Pi(x, r), and 
flavour fi(x,r), with i referring to cluster 
i. There are no constituents explicitly kept 
track of, and there is no memory referring 
to the production process. 

Instead of using a grid, we will use a dif- 
ferent but equivalent method: we assign a 
"critical volume" Vo to all particles - again 
at given time r. Domains of high energy 
density correspond to overlapping volumes 
(see fig. HI). So the task here is to find con- 
nected objects with at least pairwise over- 
lap of individual volumes (in fig. [I], we 
find four such objects). These objects are 
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Figure 2: The collision zone (in 
the t-z plane) of an AA collision. 
Each dot represents an origin of a 
string evolution. 

now considered as quark matter clusters (or 
droplets). It should be noted that the "crit- 
ical volume" Vo is not the usual volume 
47ri?!/3 of a hadron, it is rather defined by 
the requirement that if hadronic matter is 
compressed beyond Vo per hadron, the in- 
dividual hadrons cease to exist and quark 
matter is formed. So Vo is less than AnR\/?>. 
Rather than the critical volume Vo, we usu- 
ally use the critical energy density 



with m/j being (m n + m p )/2. This criti- 
cal energy density Eo (or equivalently Vo) is 
the crucial parameter of our approach. We 
will investigate properties of the system as 
a function of Eo, and based on that we will 
try to find the "realistic value" of Eo. 

For the first stage we use the indepen- 
dent string model, to be more precise the 
basic VENUS model without final state in- 
teractions QTJ] . The elementary elastic inter- 
action is Pomeron exchange, inelastic scat- 
tering and in particular particle production 
is treated via the optical theorem and by 
using the AGK cutting rules. An elemen- 
tary inelastic process is then represented by 
a "cut Pomeron", which amounts to colour 



Figure 3: Trajectories of string 
fragments. 

exchange and the creation of two strings. 
In this formalism, by using Carlo methods, 
one can determine projectile and target nu- 
cleons which interact with each other, and 
how many colour exchanges occur per NN 
interaction. In fig. |2], we show a typi- 
cal event: we consider a projection to the 
t-z plane, nucleon trajectories are repre- 
sented by straight lines, interactions are rep- 
resented by dots. Some nucleons interact 
(participants), some survive the interaction 
zone (spectators). Each dot, representing 
interaction, is a point of string formation, 
or, in other words, the origin of a string 
evolution. What happens after a string for- 
mation point? We use the standard proce- 
dure of classical relativistic string dynamics 
and decays In fig. |3] a typical example 
of the space-time evolution is shown. The 
upper rectangles represent produced parti- 
cles (hadrons and resonances), the arrows 
indicate particle trajectories. Remarkable is 
the strict ordering of the directions, being a 
consequence of the covariant string breaking 
mechanism. 

Being able to construct in a first stage, 
event by event, particle trajectories defined 
by their origins in space and the four- 
momenta, we can proceed to stage two, the 
analysis of energy densities and cluster for- 
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mation at fixed time. The question is, what 
we mean by fixed time, which frame we are 
using. This is going to be discussed in the 
following. 

Crucial for the whole approach is a corre- 
lation between the rapidity 
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and the space-time rapidity 
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with E and p z being energy and longitu- 
dinal momentum of a particle, t being the 
time, and z being the longitudinal coordi- 
nate. The variable ( is like an angle, con- 
stant ( are straight lines through the ori- 
gin in the t-z plane. In order to investigate 
a correlation between y and (, we perform 
a VENUS simulation for central S+S colli- 
sions at 200 GeV. We measure, for different 
times, the average rapidity y of produced 
particles in small cells at f± = 0, func- 
tion of z. We average over many events. We 
observe, as shown in ref. ||, to a very good 
accuracy 

V = C, (4) 

with deviations only around z ~ t, due 
to the finite size of the nuclei. We may 
parametrize our findings as 



y(t, z) 



( T for z < z T 

C(t, z) for z T < z < z p 

( p for z > z p 



(5) 



with the boundaries z p and z T given as 

z p (t) = -z T (t) = at , (6) 
with a < 1, and with 

C p = -C T = C(t,z p (t)) . (7) 

So we have an "inner region" with Bjorken- 
type behaviour y — £, where all the particle 
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Figure 4: Space-time evolu- 
tion of nucleus-nucleus scattering. 
The dotted line represents con- 
stant proper time (big dots: hy- 
perbola, small dots: tangent). 

momenta point back to the origin, and an 
"outer region" with parallel velocity vectors. 
Correspondingly, the dotted line in fig. |J a 
hyperbola in the inner region and the tan- 
gents at z p and z T , represent equal proper 
time t (on the average). In this way, ob- 
viously, also an average comoving frame is 
defined. The hypersurface defined by the 
hyperbola/tangents (dotted line) together 
with arbitrary xj_ is called r-hypersurface. 
We are now in a position to specify the 
frame for interactions: we investigate densi- 
ties (or overlap) at constant r, which means 
on r-hypersurfaces. 

Having specified the frame and corre- 
spondingly the time coordinate r, we have 
to introduce a useful longitudinal coordi- 
nate. We use the "proper length" 



dz*, 



with an integration at constant r, and with 
dz* being a longitudinal length in the av- 
erage comoving frame defined by the r- 
hypersurface. So we have 

dz' , . 

o cosh Q 

withC = C(M) = 0.5 \n(t+z)/(t-z) in the 
inner region and ( = ( T / p in the outer re- 
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gion. In the inner region we have the simple 
relation 

8 = t( (10) 

between the length s and the "angle" (. 

To specify the geometrical properties of 
particles or clusters, we use the variables r, 
s, and r_j_ = (r x ,r y ). At given r, a particle 
or cluster is considered as cylinder in s, r x , r y 
space, with the axis along the s-axis. The 
object is characterized by a lower and upper 
value of s, s\ and s 2 , and a transverse radius 
rj_. We also use 

s:= (si + s 2 )/2 (11) 

and 

As := (s 2 - Sl )/2 . (12) 

We are now in a shape to construct high 
density domains, to be considered as quark 
matter cluster. As discussed above, these 
domains are constructed, as in percolation 
models, by investigating geometrical over- 
laps of individual objects. There are two 
types of objects: initially, we have only 
hadrons and resonances, later also clusters 
contribute, which have been formed earlier. 
In any case, at a given time r, the ob- 
jects are considered cylindric in (s,r x ,r y )- 
space. Connected overlapping objects in 
this space define clusters. Such a cluster 
has in general a very irregular shape, which 
is "smoothened" in the sense that this com- 
plicated shape is replaced by a cylinder of 
the same volume, the same length s 2 — si, 
and the same center s. 

Starting at some initial time r (presently 
1 fm/c), we step through time as r i+1 = 
Ti + St, constructing clusters at each time 
Tj. For the time evolution of clusters, we 
presently assume purely longitudinal expan- 
sion, 

AC(r m ) = ACfa). (13) 

Crucial for our percolation approach is 
the initialization, i.e., the volume Vq as- 
signed to the hadrons and resonances. The 



"critical" volume Vo, or equivalently the 
"critical" energy density eo := m^/Vo, is the 
percolation parameter. To say it again, the 
volume Vo is not the usual nucleon volume 
- it is the minimum volume per particle for 
hadronic matter to exist. In order to avoid 
confusion we therefore prefer to use the crit- 
ical energy density eo as percolation parame- 
ter. There are two things we are going to do: 
first, we will investigate, for a given reac- 
tion, the dynamics SIS db function of e , then 
we will try to find a "realistic" value of eo 
(denoted as ej), by comparing with data. 

The last and most difficult topic to be 
discussed is the hadronization. The power 
of our percolation approach is that any 
hadronization scheme can be plugged into 
our approach and tested in a very detailed 
fashion. This is what we plan for the fu- 
ture. Currently, we present a very simple 
scenario, which can be implemented quite 
easily. Since our cluster expands, the energy 
density decreases, and drops at some stage 
below eo- Per definition, we hadronize the 
cluster at this point instantaneously. Since 
the clusters turn out to have essentially lon- 
gitudinal shape (As ^> Ar±), we proceed as 
follows. The cluster is cut into many short 
pieces in s, all of them having the same mass 
m with the requirement of m being around 
some parameter m seg . The small clusters 
then decay isotropically according to phase 
space ||. So this is essentially the decay of 
many fireballs at different rapidities. 

Our approach is, apart from the relativis- 
tic expansion, a typical percolation prob- 
lem. In general, one investigates the value 
for an order parameter P as a function of the 
percolation parameter. For infinite systems, 
one finds second order phase transitions, 
with P being zero in one phase and nonzero 
in the other. But also for finite systems 
a characteristic behaviour survives. In our 
case, we expect the following: for large val- 
ues of the percolation parameter eo (or small 
Vq), we expect the distribution of cluster 
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Figure 5: Distribution of cluster 
volumes for different values of the 
critical energy density (CED) e . 



sizes being peaked at small sizes, dropping 
very fast. For small eo (large Vq), it is most 
likely to find large clusters, so the distribu- 
tion will peak at large sizes. There should 
be a transition region, around some "perco- 
lation transition value" e^ T , with large fluc- 
tuations in cluster sizes. The question is 
what is the value of e^ T , in particular re- 
lated to the realistic value eg (to be deter- 
mined) . 

In fig. ||, we show the results. We plot 
the distribution of cluster volumes V for dif- 



ferent values of the critical energy density 
(CED) e - The numbers are not normal- 
ized, we show the number of clusters per 
volume bin AV^, found in 1200 simulations; 
the bin sizes for the three distributions are, 
from top to bottom: AV=10 fm 3 , 20 fm 3 , 
60 fm 3 . From the grey-scales, we can also 
read off how the clusters are distributed in 
energy density. We show only results for 
r = 2 fm/c, since volumes and energy den- 
sities turn out to scale in a trivial manner, as 
V ~ r and e ~ r _1 , so the distributions for 
different r look similar up to a scale trans- 
formation. 

We observe exactly what we expected: for 
Eo = 1 GeV/fm 3 , the distribution peaks at 
small values of V, dropping very fast with 
increasing volume V. Reducing the CED to 
0.50 GeV/fm 3 , the distribution gets wider, 
roughly by a factor of two. Although, as 
for Eq = 1 GeV/fm 3 , very small clusters 
are favoured, the fluctuations are consid- 
erably larger. Reducing Eq further to 0.15 
GeV/fm 3 (= nuclear matter density), we ob- 
serve a drastic change in the distribution, a 
maximum at large values of V emerged. It is 
possible to introduce "a maximum volume" 

V r 

given as 



max, which is for a longitudinal expansion 



Vrr, 



rAyAj 



(14) 



with the rapidity range Ay and the trans- 
verse area A±. For central S+S at 200 GeV, 
we obtain for r = 2 fm/c: 



600fm 3 



(15) 



For e = 0.15 GeV/fm 3 , the tail of the 
^-distribution almost reaches out to this 
value, and therefore this value of Eo is close 
to the value Eq T for the percolation tran- 
sition (if the probability P for V max to be 
observed is the order parameter). 

It is difficult and not yet conclusive to find 
a "realistic" value Eq for e . The best value 
at present is Eq = 1 GeV/fm 3 , which pro- 
vides good results concerning proton, pion, 
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and strange particle yields for S+S and 
S+Ag scattering. We expect further confir- 
mation (or not) from investigating particle 
correlations, because the value of e affects 
quite strongly "source radii" deduced from 
correlation functions. Corresponding work 
is under progress. 

Taking our best value, 6$ = 1 GeV/fm 3 , 
the upper plot of fig. |5| represents the real- 
istic world (the lower plots are just mathe- 
matical exercise). Although here it is most 
likely to produce just small clusters (mainly 
hadrons and resonances), there is neverthe- 
less a reasonable probability to form big 
clusters, with our statistics of 1200 simula- 
tions up to 65 fm 3 . The question is, whether 
such "miniplasmas" can be "isolated" in 
event-by-event experiments. 

To summarize, we have introduced a per- 
colation approach for the final stage (r > 1 
fm/c) of ultrarelavistic energies. This ap- 
proach allows for the formation of quark 
matter clusters in case of high energy den- 
sities. We analyzed the distribution of clus- 
ter sizes for central S+S collisions at SPS 
energies (200 GeV). Although, for a realis- 
tic value of our percolation parameter e , 
it is most likely to form only small clusters 
(mainly hadrons and resonances), there is a 
reasonable rate to form large clusters. So at 
SPS energies, it is crucial not to investigate 
average events, but to pick out events with 
large clusters, which occur due to geometri- 
cal fluctuations. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9403010vl 
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